Universitatea Tehnica a Moldovei Rezistenta Materialelor cu MATHCAD CHISINAU - 2002 Ministerul Învatamântului al Republicii Moldova Universitatea Tehnica a Moldovei Catedra Rezistenta Materialelor V BALAN REZISTENTA MATERIALELOR CU MATHCAD CHISINAU UTM 2002 Autorii acestei lucrari si-au propus sa puna la dispozitia celor interesati cunostinte minime necesare despre posibilitatile, performantele si manevrarea pachetului de programe MATHCAD si exemple de aplicare ale lui pentru rezolvarea problemelor de Rezistenta Materialelor Prezentul material didactic se adreseaza studentilor si doctoranzilor de la toate facultatile Universitatii Tehnice a Moldovei Redactor responsabil: prof univ , dr hab Vasile Marina (c) UTM, 2002 Prefata În anii 60-70 ai secolului trecut pentru rezolvarea unei probleme matematice cu ajutorul calculatorului era necesar de a scrie, verifica sintaxa si semantica, de a testa si de a pune în functie un program în unul din limbajele de programare BASIC,FORTRAN,PASCAL,C cunoasterea carora necesita un efort esential suplimentar intelectual si de durata din partea solicitantului Practic utilizatorul trebuia sa fie si programator Aceasta incomoditate frâna considerabil aplicarea calculatorului în activitatea didactica si cea inginereasca În aceste conditii, la începutul anilor 80 a devenit tot mai clar, ca produsele sofware create trebuie sa permita descrierea si rezolvarea diferitor probleme matematice fara ca utilizatorul sa posede cunostinte profunde în programare Perfectionarea metodelor numerice de calcul, aparitia, dezvoltarea si ieftinirea calculatoarelor personale, elaborarea programelor de calcul performante si simplificarea accesului la ele, interfetele prietenoase utilizatorului au condus la o explozie informationala fara precedent în istoria omenirii, la posibilitatea rezolvarii multor probleme matematice cu efort minim, care pâna nu demult pareau extrem de complicate La momentul actual, pentru calcule matematice tot mai frecvent se utilizeaza nu limbajele traditionale de programare, ci pachete de programe specializate, cum sunt MATHEMATICA, MATLAB,GAUSS,MATHCAD etc Din aceste produse cel mai popular este MATHCAD-ul, elaborat (1981) si dezvoltat de compania Math Soft Inc , Cambridge, Massachusetts (SUA) Aceasta popularitate se datoreaza în primul rând faptului, ca este unicul sistem de automatizare a calculelor matematice, în care descrierea rezolvarii problemei matematice se face asa cum o scriem pe hârtie Astfel el este accesibil si savantilor, si inginerilor, si studentilor, si elevilor Din aceasta cauza el se mai numeste "supercalculator de buzunar" Lucrarea de fata îsi propune sa ofere cititorului cunostinte initiale despre lucrul cu programul MATHCAD, ce constituie continutul primei parti, si aplicarea lui pentru rezolvarea problemelor de Rezistenta Materialelor în partea a doua a lucrarii Modul de lucru cu programul este prezentat pas cu pas, astfel ca sa fie accesibil si celor care nu au pregatire si experienta în domeniul utilizarii calculatoarelor Prima parte contine prezentarea generala a programului, trecerea în revista a comenzilor principale cuprinse în meniul principal, tipuri de date si functii, reprezentari grafice, calcule simbolice, procesorul de texte, redactarea si imprimarea rezultatelor, derivare si integrare analitica si numerica, rezolvarea ecuatiilor algebrice si diferentiale În partea a doua, pentru consolidarea cunostintelor obtinute si prezentarea practica a posibilitatilor si facilitatilor oferite de acest program se rezolva o serie de probleme traditionale ale cursului de Rezistenta materialelor, ce de obicei constituie continutul lucrarilor grafice de calcul, si probleme netraditionale, care nu se propun spre rezolvare din cauza complexitatii si volumului de calcul numeric manual excesiv de mare În aceasta lucrare se face doar o initiere în MATHCAD, pentru a arata posibilitatile de calcul si usurinta cu care el se manevreaza O descriere mai amanuntita poate fi gasita în sutele de carti dedicate acestui program Autorii 1 INITIERE ÎN MATHCAD 1 1 PREZENTAREA PROGRAMULUI Pachetul de programe MATHCAD este destinat solutionarii problemelor matematice, atât numeric, cât si simbolic si este elaborat si comercializat de compania Math Soft Inc (SUA) Acesta se poate instala pe orice calculator personal ocupând 16 MB spatiu pe disc Ruleaza sub sistemul operational Windows 95,98,2000 În raport cu alte produse similare, mentionate anterior, se remarca prin simplitate si pretentii hard reduse Posibilitatile si facilitatile oferite de ultima versiune MATHCAD 2001sunt: * Expresiile matematice se tasteaza asa cum le scriem pe hârtie: linia fractionala "-", diferentierea , integrarea , calculul limitei functiei , etc ; * Redactarea, verificarea sintaxei, compilarea si calculul se efectueaza simultan, ce esential simplifica rezolvarea problemelor De exemplu, în procesul rezolvarii problemei se poate construi graficul unei functii, analiza careia ne poate ajuta la evidentierea greselilor sau la determinarea urmatorilor pasi pentru solutionarea problemei; * Rezolvarea numerica, iar în unele cazuri si analitica a ecuatiilor si sistemelor de ecuatii algebrice liniare si neliniare (pâna la 200 ecuatii); * Rezolvarea ecuatiilor diferentiale obisnuite liniare si neliniare de diferite ordine; * Rezolvarea unor ecuatii diferentiale cu derivate partiale; * Procesorul grafic permite realizarea grafica a liniilor, suprafetelor, figurilor spatiale în coordonate carteziene, polare, cilindrice si sferice; * Existenta procesorului de texte permite crearea articolelor, rapoartelor, lucrarilor de curs gata de publicare; * Facilitati de cautare a explicatiilor prin index, descrierea amanuntita a principalelor actiuni MATHCAD, un help - context senzitiv; * Posibilitatea de a transmite o parte a problemei în mediul altui program în statica (prin memoriul tampon Clipboard) sau dinamica (OLE-tehnologii) si acolo de a o rezolva cu o eventuala revenire în mediul MATHCAD; * Pe lânga produsul MATHCAD propriu-zis sunt elaborate manuale electronice cu exemple pentru diferite disciplini: rezolvarea ecuatiilor diferentiale, statistica, termodinamica, rezistenta materialelor, etc (pot fi procurate cu plata suplimentara) În dependenta de complexitatea problemelor rezolvate se produc trei variante: - MATHCAD Standard - versiune simplificata, care se aplica în scopuri didactice; - MATHCAD Profesional - versiune profesionala, orientata pentru profesionisti; - MATHCAD Premium - versiune profesionala, destinata matematicienilor profesionisti si savantilor 1 2 ECRANUL MATHCAD Interfata utilizatorului programului MATHCAD este asemanatoare cu interfata programului WORD sub Windows, si consta din urmatoarele bare (Fig 1) enumarate de sus în jos: ==> Bara de titlu; ==> Bara de meniuri (File, Edit, View, Insert, Format, Math, Symbolics, Window, Help); ==> Bara cu instrumente Standard; ==> Bara cu instrumente de formatare Formatting ==> Bara operatiilor matematice Math (verticala, la extremitatea stânga a ecranului); ==> Fereastra de editare (aici se descrie problema); ==> Linia de stare Status Bar Fig 1 Ecranul MATHCAD Bara de titlu contine denumirea fisierului deschis La extremitatile barei sunt prezentate butoanele tipice de gestionare a ecranului (reducere, marire si închidere) La crearea unui document nou denumirea implicita este Untitled:1, care poate fi ulterior modificata Bara de meniuri si submeniurile corespunzatoare contin toate comenzile necesare pentru a efectua calcule si a edita un document Activarea comenzilor poate fi efectuata cu mouse-ul sau cu tasta Alt, tastele sageti si tasta Enter Ajustarea lor este intuitiv clara si asemanatoare cu ajustarea altor aplicatii ce activeaza sub sistemul operational Windows Comenzile care nu sunt disponibile în momentul dat sunt de culoare gri Bara cu instrumente Stadard, bara de formatare Formatting, bara operatiilor matematice Math, dubleaza comenzile sau blocurile de comenzi cele mai frecvent utilizate ale meniului principal Desenele de pe butoanele acestor bare indica operatiile si actiunile care vor fi efectuate si sunt intuitiv clare Fereastra de editare ocupa cea mai mare parte a ecranului Aici se rezolva si se redacteaza problema Initial semnul (locul) de înserare a caracterelor are forma unei cruce mici de culoare rosie, care la deschiderea documentului se afla în partea stânga-sus a ferestrei de editare (fig 1) Deplasarea semnului de înserare se efectueaza cu tastele sageti Va puteti deplasa si cu tastele Tab sau Enter Putin mai la dreapta de centrul ferestrei de editare observati o linie verticala care delimiteaza paginile desfasurate pe orizontala Daca deplasati semnul de înserare pe verticala în jos veti observa linii orizontale, care delimiteaza paginile pe verticala Deplasarile rapide în cadrul documentului se pot efectua cu ajutorul barelor de derulare situate în partea de jos si din dreapta ferestrei de editare Paginile din extrema stânga a ferestrei sunt pagini de baza ale documentului, iar celelalte se folosesc pentru comentarii sau sunt parti ale documentului care nu este necesar sa le editam De exemplu, unele calcule intermediare, care pot fi destul de voluminoase Linia de stare contine informatii despre starea actuala a programului Daca sistemul nu efectueaza nici o operatie ea afiseaza mesajul Press F1 for Help (tasteaza F1 pentru ajutor) La scrierea relatiilor matematice bara cea mai utilizata este bara operatiilor matematice Math (Fig 2) Daca lipseste pe ecran o puteti însera cu comenzile View, Toolbars, si efectuati click pe inscriptia Math Fig 2 Bara operatiilor matematice Aceasta bara este flotanta si poate fi ancorata sus, pe marginile stânga, dreapta sau la baza ecranului Pe ecranul prezentat în fig 1 ea este fixata pe marginea stânga Bara contine noua iconite, un click cu mouse-ul pe fiecare din ele deschide noua palete (fig 3): 1 Calculator Toolbar 2 Graph Toolbar 3 Vector and matrix Toolbar 4 Evalution Toolbar 5 Calculus Toolbar 6 Boolean Toolbar 7 Programming Toolbar 8 Greek Symbol Toolbar 9 Symbol keyword Toolbar Menirea lor este clara din denumire Click pe butoanele paletelor conduce la aparitia în locul semnului de înserare a cifrei, literei alfabetului grecesc sau a unui sablon, ce poate contine un operator matematic si dreptunghiuri mici pline de culoare neagra, care trebuie completate Aceste dreptunghiuri se numesc pozitii marcate ale sablonului De exemplu, click pe butonulva duce la aparitia sablonului integralei definite Daca completam pozitiile marcate si efectuam click pe butonul "=" al paletei "Calculator" sau direct pe tastatura, obtinem: Trecerea de la o pozitie marcata la alta se face cu tasta sau click stânga cu mouse-ul pe pozitia dorita Daca executati click pe integrala numai ce calculata veti observa ca expresia este încadrata într-un dreptunghi numit regiune de lucru Vizualizarea mai evidenta a regiunilor de lucru se poate efectua cu comanda View Regions a meniului Edit În acest caz fundalul ecranului devine de culoare gri, iar regiunea de lucru de culoare alba Pozitia regiunii de lucru în cadrul documentului este determinata de un punct ce se observa în partea stânga a regiunii de lucru Fig 3 Paletele matematice Pentru a muta regiunea de lucru într-o alta pozitie, fixa-ti cursorul mouse-lui pe una din liniile dreptunghiului de selectie a regiunii, indicatorul mouse-lui se va transforma într-o palma deschisa, si tinând apasat butonul stâng deplasa-ti mouse-ul în pozitia noua Daca distanta de deplasare este mare, trebuie sa va folosi-ti de comenzile Cut si Paste din bara cu instrumente Standard Paletele ocupa o parte considerabila a ecranului Din fericire, la rezolvarea unei parti a problemei nu se folosesc simultan toate paletele Le puteti închide cu click pe butonul din dreapta sus a fiecarei palete Utilizatorii avansati de obicei nu folosesc paletele în procesul de lucru Fiecare buton al paletelor are un analog selectat de la tastiera De exemplu, click pe tasta conduce la aparitia sablonului derivatei functiei Tastele pentru comenzile cele mai utilizate sunt indicate în tabela 1 Tabelul 1 Operatia Simbol pe ecran Se tasteaza Adunare Scadere Înmultire Împartire Egal (de evaluare) Egal (de definire locala) + Egal (de definire globala) + Egal (în ecuatii si relatii) + Paranteze duble Ridicare la putere + Radacina patrata Radacina de ordin superior + Diferentiere + Diferentiere de ordin superior + + Integrala nedefinita + Integrala definita + Factorial + Valoare absoluta + Suma dupa un indice + Suma pe domeniu + + Produs dupa un indice + Produs pe domeniu + + Mai mare decât > Mai mic decât Mai mare sau egal + Mai mic sau egal + Diferit + Indice inferior (textual) (exemple) Indice inferior (matematic) (exemple) Indice superior (extragerea unei coloane dintr-o matrice) (exemplu) + Indice superior textual (indice prim) (exemplu) + Determinantul matricei Transponarea vectorului sau matricei + Vectorizare (exemplu) + Transportarea unei parti a expresiei pe urmatorul rând (exemplu) + Crearea regiunii de text Calculele în MATHCAD se executa de la stânga la dreapta si de sus în jos Astfel, pozitiile regiunilor de lucru în cadrul documentului determina consecutivitatea efectuarii calculelor 1 Selectam MA=q*L2 si afisam rezultatul 1 3 OPERATORI DE DEFINIRE SI OPERATORI DE CALCUL Semnul " = " scris pe hârtie, în dependenta de context, are sensuri diferite: el poate defini o functie, poate atribui unei variabile sau parametru o valoare sau câteva valori (vector sau matrice), poate afisa valoarea unei variabile, functii, vector sau matrice, poate însemna efectuarea evaluarii unei expresii, sau egalul într-o ecuatie Evident ca programul MATHCAD trebuie sa înteleaga ce se are în vedere Pentru definiri si atribuiri se utilizeaza operatorul ":=", care se obtine pe ecran cu tasta , pentru afisari si evaluari operatorul " = ", se obtine cu tasta , pentru ecuatii operatorul " = ", se obtine cu combinatia de taste + Acesti operatori se pot afisa si de pe paletele "Calculator" si "Boolean" (fig 3) Observati ca semnul egal pentru ecuatie este mai accentuat decât semnul egal pentru afisare Dupa cum s-a mentionat calculele în MATHCAD se efectueaza de la stânga la dreapta si de sus în jos Deci, definirea datelor initiale trebuie sa se faca înainte de a scrie relatiile si ecuatiile de calcul Dar în unele cazuri este comod de a scrie datele initiale dupa relatiile de calcul MATHCAD-ul poseda un asa operator "=", numit operator de definire globala Se obtine tastând + , sau de pe paleta "Evaluation" Definirea astfel efectuata este valabila pe întreg documentul Cu alte cuvinte, calculul se poate efectua înainte de definire Aplicarea acestor operatori va deveni mai clara din exemplele prezentate ulterior 1 4 CALCULUL EXPRESIILOR ARITMETICE SI ALGEBRICE Asa dar, în fata avem ecranul MATHCAD cu denumirea implicita a documentului "Untitled 1" Pentru început vom da denumire primului nostru document si îl vom salva Pentru asta efectuam urmatorii pasi: deplasam cursorul mouse-lui pe bara de meniuri, efectuam click pe meniul File, alegem optiunea Save; în fereastra aparuta Save alegem dosarul dorit; în câmpul File Name introducem numele fisierului, sa zicem "Ex1", si efectuam click pe butonul Save Documentul a obtinut denumire si este salvat Este bine sa alegeti formatul paginii si sa stabiliti marginile ei Pentru asta efectuati comanda File, Page Setup Pe ecran va aparea fereastra de dialog Page Setup (Fig 4) În caseta derulanta a câmpului Size alegeti formatul, sa zicem A4 si stabiliti marginile paginii în câmpurile Top(sus), Bottom(jos), Left(stânga), Right(dreapta) Apasati OK Implicit semnul de înserare " + " de culoare rosie se gaseste în partea stânga-sus al primei pagini de lucru Deplasati-l putin mai jos si spre dreapta cu tastele sageti Pentru a vedea mai bine regiunile de lucru executati click pe meniul View si alegeti optiunea Regions Ecranul a devenit de culoare gri Acum suntem pregatiti pentru a efectua primele calcule Fig 4 Formatarea paginii Exemplu 1 Vom calcula valoarea expresiei aritmetice Pentru aceasta efectuati urmatoarele actiuni: 1 Selectati consecutiv de la tastiera , , Odata cu introducerea cifrei " 2 " s-a deschis regiunea matematica de lucru, iar semnul de înserare s-a transformat într-un cornier de culoare albastra orientat spre stânga, care cuprinde cifra Dupa ce ati selectat punctul zecimal si cifra " 1 " observati ca cornierul s-a extins pe tot numarul, astfel de parca îl sustine Cu tasta se sterg semnele din interiorul cornierului, iar cu tasta din exteriorul lui Schimbarea orientarii cornierului se poate face cu tasta 2 Selectati operatorul de la tastiera sau de pe paleta "Calculator" Ati obtinut pe ecran: 3 Apasati consecutiv , + , Acum cursorul se afla pe pozitia exponentei Pentru a coborî de pe aceasta pozitie apasati tasta de spatiu Cursorul s-a extins pe termenul 32 4 Selectati , , Apasati de trei ori Cursorul s-a extins pe toata expresia 5 Selectati , , + , , + Pe ecran a aparut sablonul radacinii de ordin superior Apasati consecutiv , , 6 Extindeti cursorul pe tot numitorul prin apasarea tastei de doua ori Selectati Numitorul s-a încadrat în paranteze Apasati Cursorul a cuprins si parantezele Selectati , 7 Declansam procesul de calcul cu tasta , sau cu butonul "=" de pe paleta "Calculator" Am obtinut: 8 Apasati pentru a parasi regiunea matematica de lucru Semnul de înserare s-a deplasat în jos 9 Tuturor numerelor în MATHCAD le este asociat implicit stilul de formatare "Constants", caruia la rândul sau îi este asociat fontul "Times New Roman" si dimensiunea "10" Le vedeti în bara "Formatting", în partea de sus a ecranului 10 Acest stil poate fi modificat de utilizator Efectuati click-stânga pe sageata de derulare a dimensiunii fontului "Font Size" si din lista aparuta selectati dimensiunea "12" Daca activati unul din butoanele "Aldin", "Cursiv" sau/si "Underline" puteti obtine numere cu aldine, cursiv sau/si subliniate Mentionam, ca stilul "Constants" stabilit se va aplica tuturor numerelor, rezultatelor numerice si chiar exponentelor numerice a documentului curent 11 Frecvent se întâmpla ca doriti sa reamplasati într-o pozitie noua regiunea de lucru Efectuati click pe regiunea de lucru Stabiliti cursorul mouse-lui pe una din laturile dreptunghiului ce încadreaza regiunea Cursorul se transforma într-un semn asemanator cu o palma deschisa Apasati butonul stâng al mouse-lui si tinându-l apasat deplasati-l în pozitia dorita Eliberati butonul 12 Pentru o vizualizare mai evidenta folositi-va de butonul "Zoom" (Scara) Alegeti din lista derulanta scara "150%" Exemplu 2 Calculam expresia , unde a=5 si ß=30º 1 Dupa cum am mentionat calculele în MATHCAD se fac de la stânga la dreapta si de sus în jos Deci, pentru a determina valoarea marimii X este necesar în prealabil de atribuit parametrilor "a" si "ß" valori numerice Aplicam operatorii de definire locala Stabilim cu mouse-ul sau cu tastele sageti pozitia semnului de înserare, de unde vom începe a efectua calculele Apasam consecutiv , + , , Din paleta alfabetului grec "Greek" apasam "ß", apoi + , , , , , , Am obtinut: 2 Atribuim marimii "X" valoarea expresiei din partea dreapta a semnului "egal" Apasam + , + , si din paleta "Calculus" butonul sumei 3 Selectati , , , + , , , din paleta "Greek" butonul "p", , , , + , , , , , , , Programul recunoaste unele din constantele matematice inclusiv "p" si "e" Puteti sa va convingeti de acest lucru, daca afisati valorile lor Pentru toate functiile argumentul trebuie încadrat în paranteze 4 Apasati , , + , , , Pentru a selecta termenul doi este necesar de a extinde cornierul de selectie pe tot primul termen Apasati de trei ori tasta Dupa aceasta introducem operatorul "+", iar de pe paleta "Calculator" apasam butonul logaritmului natural "ln" 5 Introducem , + , Pentru a lua cubul logaritmului natural extindem cornierul de selectie pe tot logaritmul apasând de doua ori tasta Selectam + , Apasam , , butonul tangentei "tan" de paleta "Calculator", litera "ß" de pe paleta "Greek" si 6 Afisam valoarea marimii X Pentru aceasta stabilim o pozitie noua a semnului de înserare si selectam + , Am obtinut: 7 Toate calculele programul MATHCAD le face cu exactitatea de 16 cifre semnificative Implicit formatul de afisare pe ecran si de editare a rezultatelor numerice prevede trei cifre zecimale Dar el poate fi foarte simplu modificat Efectuati dublu-click cu butonul stâng al mouse-lui pe orice rezultat numeric Imediat apare fereastra Result Format In câmpul "Number of decimal places" (numarul de cifre zecimale) stabiliti exactitatea de afisare dorita Pentru marimea X am stabilit exactitatea maxima de 15 cifre zecimale Observati ca programul propune câteva formate ,care sunt indicate în caseta de dialog "Format" Cel mai uzual este formatul "General" Aceasta fereastra propune înca câteva optiuni Încercati-le Modificarea formatului numeric poate fi realizat nu numai local (pentru un singur rezultat), dar si global (pentru toate rezultatele numerice din document) Efectuati click undeva în afara regiunilor de lucru Selectati comanda Rezult din meniul Format Pe ecran va aparea aceeasi fereastra Result Format, dar modificarile stabilite în ea vor fi valabile pentru tot documentul 8 Mai jos de definirea marimii X afisati înca odata valoarea ei, si deplasati-o în sus Programul printr-un mesaj de eroare ne anunta, ca aceasta variabila sau functie nu a fost prealabil definita 9 Efectuati click pe una din literele expresiei matematice Pentru toate literele din expresiile matematice programul propune implicit stilul Variables cu fontul Times New Roman de dimensiunea 10 din bara de formatare Formatting Dimensiunile numerelor si literelor este de dorit sa fie aceleasi Schimbati dimensiunea fontului din 10 în 12 si stabiliti stilul accentuat (Bold) al literelor Am obtinut: Exemplu 3 Valorile numerice în calculele ingineresti sunt urmate de unitati de masura MATHCAD-ul permite în unele cazuri lucrul cu unitatile de masura Unitatile de masura se aleg si pot fi vizualizate din fereastra de dialog Insert Unit (Fig 5), care se afiseaza cu comanda Unit a meniului Insert Ne vom familiariza cu aplicarea lor 2 Selectati consecutiv , + , , , + , , , , , , , , , + , , , , Ati observat ca pentru a atasa unui parametru sau variabile unitati de masura se tasteaza dupa valoare Am obtinut: 3 Unitatile de masura pot fi înserate nu numai de la tastatura, dar si din fereastra Insert Unit î Fig 5 Fereastra unitatilor de masura 4 Se pot defini si alte unitati decât cele oferite de listele din Insert Unit Unitatea noua definita printr-o operatie de atribuire va fi valabila în continuare în document 5 Introduceti , , + , , , + , , , , Selectam , , , + , , , , Afisam valoarea marimii VA , , , Obtinem: 6 În ultima regiune de lucru stabiliti cornierul de selectie pe pozitia marcata din dreapta (cu tastele sageti sau click cu butonul stâng al mouse-lui) Selectati , , Rezultatul va fi afisat în kN Daca nu va sunt pe plac zerourile care apar în rezultat, chemati fereastra Result Format si deselectati optiunea Show treiling zeros 7 Rezultatul l-am obtinut în unitati de energie J(Joules) Reamintim ca 1J=1 N*m Conversiunea la N*m se face simplu Fixati cornierul pe pozitia marcata si selectati , , , Transformati rezultatul în kN*m 8 Efectuati urmatoarele calcule: 9 Programul MATHCAD are si un procesor textual, care permite introducerea textelor Amplasati semnul de înserare undeva între primul si al doilea rând a ultimului calcul Apasati de câteva ori tasta Toata partea documentului, ce se afla mai jos de semnul de înserare se va deplasa în jos Pentru a lichida spatiile libere între rânduri folositi-va de tasta Efectuati click pe la mijlocul spatiului liber obtinut Apasati optiunea Text Region al meniului Insert Pe ecran va aparea regiunea de text, iar semnul de înserare se va transforma într-o linie verticala de culoare rosie Procesorul de texte poate fi chemat si cu tasta Selectam textul Implicit pentru regiunile de text se stabileste stilul "Normal" caruia îi este atasat fontul Arial cu dimensiunea 10 Modificati-l dupa dorinta Cu parere de rau, procesorul nu permite selectarea literelor specifice alfabetului român Pentru a parasi regiunea de text efectuati click în afara regiunii Cu tasta se trece la un rând nou în regiunea de text actuala Daca doriti sa stergeti simultan câteva regiuni de lucru, imaginati-va un dreptunghi în care se încadreaza aceste regiuni Apasati butonul stâng al mouse-lui în coltul stânga-sus al dreptunghiului imaginat, tinându-l apasat, deplasati-l în coltul drept-jos Eliberati butonul Toate regiunile selectate se vor încadra în dreptunghiuri cu linii întrerupte Acum, daca amplasa-ti sageata mouse-lui în interiorul uneia din aceste regiuni, ea se transforma într-o palma deschisa Apasa-ti butonul stâng al mouse-lui si tinându-l apasat deplasati regiunile selectate într-o pozitie noua Iar stergerea se face foarte simplu: apasati tasta 1 5 VARIABILE SIR SI VARIABILE INDICE Un sir ordonat de numere în forma de progresie aritmetica identificat printr-un nume se numeste variabila sir, care se defineste cu sintaxa: Variabila sir:=val initiala,val urmatoare val finala Ratia progresiei este egala cu diferenta dintre valoarea urmatoare si valoarea initiala Exemplu Selectati , + , , , , , , , , Variabila sir este definita (fig 6) Ati observat ca pentru a obtine operatorul sirului " " se tasteaza Afisati valorile variabilei sir cu tastele , (fig 6) Variabila sir care consta din numere naturale consecutive (ratia este 1 sau -1) se numeste variabila indice (range variable) si se defineste cu sintaxa: Variabila indice:=val initiala val finala Fig 6 Exemple de variabile sir si variabile indexate Exemplu Selectati , operatorul de definire + , valoarea initiala , operatorul sirului " "cu tasta sau cu butonul "m n" de pe paleta Matrix, valoarea finala si parasiti regiunea de lucru cu tasta Afisati valorile variabilei indice cu tastele si (fig 6) Sa admitem ca aveti nevoie de a defini un sir de numere sau expresii arbitrare Cu ajutorul variabilei indice se poate defini o astfel de variabila denumita variabila indexata cu urmatoarea sintaxa: Variabilavariab indice:=expresie 1,expresie 2, expresie N sau pentru fiecare valoare a indicelui aparte: variabila valoare variabila indice:=expresie Pozitia indicelui variabilei indexate se obtine cu tasta Între valorile sirului se tasteaza Exemplu (fig 6) Selectati numele variabilei , coborâti pe pozitia indicelui inferior cu tasta , selectati numele variabilei indice care este deja definita, operatorul de definire + , introduceti numerele sau expresiile dorite separându-le prin virgule Apasati Variabila indexata este definita Între variabilele sir si variabilele indexate exista o deosebire calitativa Variabila indexata are o singura valoare pentru o valoare concreta a indicelui Apasati , , , Ati obtinut "z1=1 1" Valoare a variabilei sir este tot sirul de numere sau expresii si de aceea nu se poate accesa o valoare singulara a ei Nu se poate scrie "M:=i-1", deoarece "M" este o variabila obisnuita, iar "i" este o variabila indice Totodata, daca vom defini o variabila indexata "Mi" cu expresia "i2+1", programul nu va avea obiectii Observati ca avem acces la fiecare valoare a variabilei indexate Valorile unei variabile indexate sunt întocmai componentele vectorului cu acelasi nume Selectati si operatorul de evaluare Veti obtine sase numere numite componente ale vectorului, aranjate pe verticala si încadrate în paranteze rotungite În MATHCAD forma desfasurata a vectorilor si matricelor se prezinta în paranteze rotungite Sase si nu cinci componente ale vectorului se datoreaza faptului, ca implicit toate variabilele indice se inizializeaza cu valoarea "0", indiferent de valoarea initiala introdusa Deci implicit prima componenta a oricarui vector are indicele "0" si daca aceasta valoare (si urmatoarele) nu sunt definite în variabila indice, componentele corespunzatoare ale vectorului vor fi afisate cu valorile "0" (vezi si fig 6) Originea indicilor poate fi schimbata cu constanta programului ORIGIN, care implicit are valoarea "0" Atribuim acestei constante valoarea initiala "1" a variabilei indice "i" Obtinem: Prezenta unei variabile indexate într-o expresie presupune evaluarea acelei expresii pentru fiecare valoare din sir în cadrul unui proces iterativ 1 6 FUNCTII SI REPREZENTARI GRAFICE În MATHCAD functiile se definesc cu urmatoarea sintaxa: nume functie(arg1,arg2, ):=expresie Argumentii sunt parametri formali, care pot fi variabile numerice, vectori, matrice sau alte functii Tipul expresiei defineste tipul valorilor functiei: numeric, vectorial sau matriceal Calculul valorii unei functii într-un punct se face înlocuind parametrii formali cu valori actuale si efectuând evaluarea Daca argumentele sunt definite cu variabile sir sau variabile indexate, evaluarea functiei va furniza sirul valorilor corespunzatoare functiei Deosebesc functii standard în MATHCAD (aproximativ 250) si functii definite de utilizator De obicei numele (identificatorul) functiei se selecteaza de la tastatura, dar daca l-ati uitat sau nu-l cunoasteti va puteti folosi de fereastra de dialog Insert Function,care se afiseaza cu comanda Function al meniului principal Insert Alegeti din lista Function Category tipul functiei, din lista Function Name numele functiei si apasati butonul Insert În fereastra de editare, în locul semnului de înserare va aparea sablonul functiei selectate Completati pozitiile marcate si afisati rezultatul cu tasta Apasati butonul semnului de întrebare în fereastra Insert Function si veti obtine informatii despre functia curenta (în limba engleza) Cele mai uzuale categorii de functii standard sunt: * Functii trigonometrice, directe si inverse (Trigonometric) * Functii hiperbolice (Hyperbolic) * Functii logaritmice si exponentiale (Log and Exponential) * Transformarile Fourier (Fourier Transform) * Functii de rezolvare a ecuatiilor si sistemelor de ecuatii algebrice liniare si neliniare (Solving) * Functii de rezolvare a ecuatiilor diferentiale (Differential Ecuation Solving) * Interpolarea si extrapolarea functiilor (Interpolation and Prediction) * Functii de regresie (Regression and Smoothing) * Functii de lucru cu vectori si matrice (Vector and Matrix) * Functii discontinue (Piecewise Continuous) Exemple de functii: Reprezentarile grafice ale functiilor se obtin, daca apasati unul din butoanele paletei "Graph" (fig 3) Puteti realiza grafice plane (rectangulare si polare) si tridimensionale (în coordonate carteziene, cilindrice si sferice) Fixati punctul de înserare într-un loc comod al ferestrei de editare, având în vedere ca graficul va fi afisat în dreapta si în jos Apasati butonul "X-Y Plot" al paletei "Graph" Ati obtinut macheta graficului plan în coordonate carteziene rectangulare Completati pozitia marcata de la mijlocul axei orizontale cu variabila independenta "x" Cu tasta treceti pe pozitia marcata de la mijlocul axei verticale si selectati de pe tastatura numele functiei cu indicarea argumentului, de exemplu sin(x) Apasati pentru a parasi regiunea de lucru Ati obtinut graficul acestei functii Daca în prealabil n-ati precizat intervalul argumentului programul îl stabileste implicit între "-10" si "+10" Dar puteti sa-l modificati Efectuati click în regiunea graficului Valorile maxime si minime ale argumentului si functiei, care pot fi modificate sunt încadrate în corniere mici Observati ca poate fi modificat si domeniul variabilei independente si domeniul valorilor functiei Pozitionati cornierul de înserare cu un click sau cu tastele sageti pe aceste valori limite si schimbati-le Dreptunghiul exterior al regiunii grafice contine câteva pozitii marcate, cu ajutorul carora puteti schimba dimensiunile ei Pe una si aceeasi regiune grafica programul permite trasarea graficelor câtorva functii Pentru aceasta, dupa denumirea primei functii introduceti virgula si selectati urmatoarea functie, de exemplu cos(x) Am obtinut desenul ce urmeaza (regiunea este activata si sunt vizibile marcajele de redimensionare si limitele modificate ale argumentului si functiei) Daca parasiti zona de lucru marcajele si limitele nu vor fi vizibile Traseul functiei cos(x) este prezentat cu linie întrerupta de culoare albastra Este posibil ca formatul graficelor nu va convine Efectuati click cu butonul dreapta al mouse-lui pe regiunea grafica, din meniul contextual aparut executati click-stânga pe optiunea Format Pe ecran se va deschide fereastra Formatting Currently Selected X-Y Plot, care contine patru sectiuni Implicit se deschide prima sectiune X-Y Axes Optiunile ce apar în aceasta sectiune au urmatoarele semificatii: Log Scale - axe logaritmice (limitele axelor trebuie sa fie pozitive); Grid Lines - permite trasarea unei retele de linii verticale si orizontale; Numbered - apar etichete în dreptul marcajelor axelor de coordonate; Autoscale - autoscalarea limitelor axelor; Show Markers - permite obtinerea a doua perechi de linii verticale si orizontale pentru evidentierea unor puncte specifice ale graficului; Auto Grid - autoretea (implicit este selectata) Activati optiunea Grid Lines Dezactivati optiunea Autoscale În câmpurile Number of Grids introduceti, de exemplu cifra "4" pentru axa "x" si "6" pentru axa "y" Din Axes Style bifati optiunea Grossed Apasati butonul "?????????" (daca este instalata versiunea rusa Windows) Am obtinut formatul din fig 7 Fig 7 Afisati sectiunea Traces (fig 8) Aceasta sectiune permite formatarea fiecarei functii separat, cu selectarea liniilor Trace 1(prima functie), Trace 2(functia 2), printr-un click pe traseul corespunzator si stabilirea optiunilor dorite din listele derulante ce se obtin cu click pe triunghiurile pline a partii de jos a tabelului Pentru trasarea graficului, programul prezinta implicit argumentul în forma de variabila sir cu ratia foarte mica Ulterior se calcula valorile functiei pentru fiecare valoare a variabilei sir Între doua valori vecine a functiei se efectuiaza interpolare liniara Acest lucru poate fi observat daca vom crea variabila sir în mod explicit În acest caz facilitatile oferite de sectiunea Traces vor deveni mai clare Putin mai sus de regiunea grafica declarati argumentul în forma de variabila sir în intervalul (-2p, 2p) cu pasul de calcul 0,4p Fig 8 Deschideti fereastra de formatare si în ea sectiunea Traces (fig 8) Fig 9 Pentru trace 1 stabiliti Symbol(o's), Line(solid), Color(blk), Type(solid) si grosimea liniei Weight(2), iar pentru trace 2 - Type(bar), Color(blk) Apasati butonul "OK" Afisati valorile functiei sin(x) Ati obtinut graficul din fig 9 Daca alegeti pasul (ratia) de calcul a argumentului, destul de mic, cu optiunea bar se pot hasura diagramele eforturilor interioare în cursul de Rezistenta Materialelor (fig 10) Fig 10 Activati graficul cu un click al mouse-lui Efectuati dublu click pe grafic Va fi afisata fereastra de formatare Executati click pe butonul sectiunii Labels În câmpurile Title si Axis Labels selectati titluri pentru grafic si pentru axe Precizati pozitia titlului graficului: Above (deasupra) sau Below (sub grafic) Definitivati operatia prin selectarea butonului "OK" Activati graficul Stabiliti cursorul mouse-lui în regiunea graficului Efectuati click dreapta În meniul aparut, selectati optiunea Zoom Cu aceasta puteti obtine detalierea unei portiuni de grafic Informatii suplimentare despre trasarea graficelor functiilor le ofera meniul Help al programului În mecanica frecvent se aplica functiile standard din categoria functiilor discontinue (Piecewise Continuous) Le vom defini: F(x) - functia Heviside, (functie treapta unitara sau functie întrerupator) care are valoarea 1, daca x=0 si valoarea 0, daca x * Calculati derivatele functiei pentru unele valori ale argumentului Cu ajutorul procesorului simbolic calculul derivatelor (si a integralelor) se poate efectua si analitic Accesul la acest procesor se face cu paleta "Symbolic" sau cu comenzile meniului "Symbolics" Prezentam câteva exemple de calcul simbolic a derivatelor Consecutivitatea actiunilor: 1 Apasati butonul derivatei de ordinul unu a paletei "Calculus"; 2 Completati pozitiile marcate; 3 Apasati butonul sageata "Symbolic Evaluation" al paletei "Symbolic"; 4 Apasati tasta Imediat veti obtine functia derivata Este regretabil faptul ca procesorul numeric si procesorul simbolic actioneaza independent si astfel rezultatele calculelor simbolice nu pot fi evaluate numeric Este posibil de a calcula si derivate de ordin superior Aplicati butonul respectiv al paletei "Calculus" Se pot calcula si derivate partiale Pentru a aplica semnul derivatei partiale, efectuati click-dreapta pe functia care se deriveaza si din meniul exploziv alegeti consecutiv comanda View Derivative As, Partial Derivative Integrarea se efectueaza în acelasi mod ca si diferentierea Fixati semnul de înserare în locul unde doriti sa fie calculata integrala, apasati butonul integralei nedefinite sau a integralei definite de pe paleta "Calculus", completati pozitiile marcate a machetelor integralelor, apasati butonul-sageata a paletei "Symbolic" pentru integrala nedefinita sau înserati operatorul de evaluare pentru integrarea definita Imediat va fi afisat rezultatul calculului Limitele de integrare pentru calculul integralei definite pot fi nu numai numere dar si variabile si chiar functii Este posibil de a calcula integrale duble si integrale triple Machetele lor se obtin simplu: apasati de doua ori respectiv de trei ori butoanele integralelor de pe paleta "Calculus" 1 8 VECTORI SI MATRICE O totalitate de numere, simboluri, functii ordonate se numeste matrice De obicei matricele constau din linii si coloane Matricea cu o singura coloana se numeste vector Componentele matricei se noteaza cu indici Primul indice arata linia, iar al doilea coloana la intersectia carora se afla elementul respectiv Definirea unei matrice se face în felul urmator: * Introduceti numele matricei urmata de operatorul de definire pe care îl obtineti cu tasta ; * Apasati butonul matricei din paleta "Matrix" Pe ecran va fi afisata fereastra "Matrix"; * Introduceti numarul de linii (Rows) si de coloane (Columns) necesare, * Apasati "Insert" Va fi afisata macheta matricei Completati pozitiile marcate Treceti de la o pozitie la alta cu tasta sau efectuati click pe pozitia dorita cu mouse-ul Matricea este definita Accesul la elementele matricei se face în felul urmator: * Selectati numele matricei; * Apasati tasta Cursorul s-a deplasat pe pozitia indicelui * Introduceti numarul liniei(prima linie are numarul 0!) în care se afla elementul cautat, virgula si numarul coloanei(prima coloana are numarul 0!) Pentru a începe numerotarea liniilor si coloanelor cu alt numar, schimbati valoarea constantei programului ORIGIN cu valoarea dorita!; * Introduceti operatorul de evaluare cu tasta Valoarea componentei va fi afisata pe ecran; Vectorii se definesc în acelasi mod Unica deosebire este ca în fereastra "Matrix" numarul de coloane trebuie sa fie 1 Daca numarul de linii si coloane este mai mare ca 15 fereastra "Matrix" nu mai poate fi folosita În acest caz este necesar de a crea matricea prin intermediul variabilelor sir Prin doua variabile sir se anunta dimensiunile liniilor si coloanelor Tuturor componentelor matricei i se atribuie valori zero(ai,j=0) Componentelor nenule li se da valori Se afiseaza matricea (daca este necesar) Implicit se afiseaza doar primele 15 linii si 15 coloane Cu barele de derulare de jos si din dreapta puteti afisa orice componente Mai mult ca atât puteti afisa matricea în întregime cu ajutorul marcajelor ce înconjoara matricea Daca efectuati un click-dreapta pe matrice, selectati comanda Properties si în fereastra aparuta apasati butonul Font puteti mari sau micsora fontul valorilor componentelor Pentru formatarea rezultatelor afisati fereastra Rezult Format, apelând comanda Result a meniului Format Cu matricele în MATHCAD se pot efectua practic toate operatiile posibile, de la adunare, scadere, înmultire, calculul determinantilor, matricelor inverse si transponarea lor, pâna la calculul valorilor proprii si vectorilor proprii Pentru calculul valorilor proprii si a vectorilor proprii ale unei matrice se folosesc o serie de functii predefinite eigenvals(T) - returneaza un vector ale carui componente sunt valorile proprii ale matricei Anxn; eigenvecs(T) - returneaza o matrice ce contine vectorii proprii normalizati la lungime unitara corespunzatori valorilor proprii ale matricei Anxn Coloana "n" a matricei returnate este vectorul propriu corespunzator valorii proprii "n" obtinute cu functia eigenvals eigenvec(T,t) - returneaza un vector propriu asociat cu valoarea proprie t a matricei T Vectorul propriu este normalizat la lungime unitara Va reamintim ca în MATHCAD se aplica doua tipuri de indici inferiori: indici-text care se includ în document cu aplicarea tastei si indici matematici care se includ în document cu tasta sau cu butonul "xn" al paletei "Matrix" Valoarea indicelui matematic evidentiaza componenta corespunzatoare a vectorului sau a matricei În documentul recent prezentat indicii inferiori 1,2,3 sunt indici matematici Pot fi efectuate si operatii simbolice cu matrice 1 9 CALCULE SIMBOLICE În MATHCAD este încorporat un procesor destinat calculelor simbolice Spre deosebire de calculele numerice efectuate cu o oarecare eroare, evaluarea expresiilor cu procesorul simbolic este exacta Calculele simbolice se pot face în doua moduri: 1 cu instrumentele paletei Symbolic; 2 cu comenzile meniului Symbolics Se foloseste în majoritatea cazurilor paleta Symbolic, dar exista situatii, când este mai comod de aplicat comenzile meniului Symbolics Pentru început vom studia aplicarea instrumentelor paletei Symbolic Evaluarea simbolica cea mai simpla se face în felul urmator: Scrieti expresia ce va fi evaluata simbolic; apasati butonul evaluarii simbolice (sageata) de pe paleta Symbolic; apasati tasta Imediat obtineti rezultatul Daca procesorul nu poate evalua expresia ea va fi returnata în forma sa initiala În acest caz trebuie sa aplicati comenzile paletei Daca termenii expresiei contin numai fractii ordinare si numere întregi rezultatul va fi fractie ordinara; daca expresia va contine si termeni în forma zecimala rezultatul va fi fractie zecimala cu 20 de cifre semnificative Pentru a micsora numarul de cifre semnificative procedati astfel: * scrieti sau efectuati click pe expresie; * apasati butonul float al paletei Symbolic; * scrieti în locul pozitiei marcate numarul dorit de cifre; * apasati tasta Pentru a descompune în factori(factor), desfasura(expand), simplifica expresiile algebrice procedati în acelasi mod Singura deosebire este aceea, ca în majoritatea cazurilor este necesar de a sterge pozitia marcata si virgula ce apare dupa apasarea butonului comenzii În unele cazuri este necesar de a efectua câteva operatii asupra unei expresii Executati-le consecutiv una dupa alta Rezultatele se vor desfasura într-o singura linie, care poate fi destul de lunga Pentru a documenta lucrarea este comod de a realiza consecutivitatea de comenzi într-un singur corp Procedati în felul urmator: * scrieti expresia sau efectuati click pe ea; * apasati butonul primei comenzi de pe paleta Symbolic; * îndeplini-ti consecutiv celelalte comenzi; * sterge-ti (daca este necesar) pozitiile marcate si virgulele; * apasati tasta Evaluarea unei expresii ce contine parametri se poate efectua dupa cum urmeaza: Substituirea parametrilor în expresii se poate face si mai elegant: Cu comanda Solve a paletei Symbolics se pot obtine solutii analitice a unor ecuatii algebrice si chiar diferentiale Pentru ecuatii algebrice de grad mai mare ca 2 solutiile sunt laborioase si de obicei nu pot fi folosite în calculele ingineresti Pentru calcule simbolice intermediare care nu vor fi documentate se pot aplica comenzile meniului Symbolic Procedati în felul urmator: 1 Scrieti expresia care va fi evaluata; 2 Cuprindeti cu cursorul de selectie (cornierul albastru) toata expresia, o parte a ei sau numai variabila în raport cu care se face operatia; 3 Alegeti din meniul Symbolic comanda care doriti s-o efectuati; Implicit rezultatul este afisat putin mai jos de expresie Efectuati urmatorul exemplu de simplificare a unei expresii algebrice La dreapta sunt scrise comenzile care sunt îndeplinite asupra expresiei din stânga: Rezultatul evaluarii poate fi afisat si la dreapta expresiei initiale cu sau fara comentarii Pentru asta afisati fereastra de dialog Evaluation Style cu comanda Evaluation Style a meniului Symbolic Stabiliti optiunile Horizontally (Orizontal) si Show Comments (Afisarea comentariilor) Urmatoarele evaluari vor fi executate în corespundere cu optiunile stabilite (Fig 13) Comentariile indica denumirea operatiei efectuate Optiunea Evaluate in Place va permite sa faceti evaluarea în acelasi loc unde a fost expresia initiala Expresia initiala dispare de pe ecran Fig 13 Cu optiunile comenzii Variable a meniulul Symbolic puteti efectua diferentierea (Diferentiate), integrarea (Integrate), substituirea (Substitute), rezolvarea ecuatiilor (Solve), descompunerea în serii (Expand to series) Pentru a realiza aceste operatii este necesar în prealabil sa evidentiati variabila în raport cu care ea se va face În caz contrar operatiile solicitate nu vor fi accesibile 1 10 ECUATII ALGEBRICE NELINIARE Orice ecuatie cu o singura necunoscuta poate fi scrisa în felul urmator: f(x) = 0 Radacina ecuatiei este valoarea variabilei x pentru care expresia f(x) se anuleaza A rezolva o ecuatie înseamna a determina numarul de radacini si valorile acestora În MATHCAD, si în alte soft-uri destinate calculelor matematice, ecuatiile se rezolva cu metode iterative (apropierea de radacina ecuatiei prin aproximatii succesive), începând cu o valoare oarecare a variabilei, stabilita de utilizator Metodele iterative, numite si metode numerice de calcul, obtin solutia (solutiile) ecuatiei cu o oarecare aproximatie, si nu exact, cum sunteti obisnuiti s-o faceti În MATHCAD aproximatia de calcul este determinata de constantele TOL si CTOL al sistemului Implicit valorile acestor constante sunt 0,001 Modificarea lor se face din fereastra de dialog Math Options, care se afiseaza cu comanda Options a meniului Math Daca radacina calculata a unei ecuatii este de exemplu x=10, aceasta de loc nu înseamna ca radacina exacta este în limitele 9,999 se include cu butonul "Equal" al paletei "Boolean" sau cu combinatia de taste + 4 Se calculeaza vectorul rezultatului (radacinii) sistemului cu functia Find 5 Se afiseaza vectorul rezultatului si componentele lui, daca este necesar 6 Daca nu ati obtinut solutia necesara, sau pentru a obtine alte solutii, schimbati valorile initiale ale problemei Raspunsul la întrebarea câte solutii are sistemul de ecuatii neliniare este o problema destul de dificila si poate fi rezolvata doar prin analiza calitativa a sistemului sau prin scanarea variabilelor initiale pe tot domeniul de variatie ale lor Cresterea preciziei calculelor se poate obtine prin micsorarea valorii constantei TOL si CTOL Ecuatiile pot fi scrise si asa: În acest exemplu, sistemul de ecuatii are doua solutii, iar în cel precedent sase! Încercati sa le determinati Cu procedura Given-Find se pot rezolva si sisteme de ecuatii ce contin restrictii de variatie a variabilelor La rezolvarea ecuatiilor neliniare deseori apar situatii, când procedurile root si given-find nu gasesc sau nu sunt în stare sa gaseasca solutii De exemplu, ecuatia din fig 15 are doua radacini, dar daca stabilim valoarea initiala a variabilei x = -1, procedura given-find nu gaseste nici una Rezultatul se afiseaza cu culoare rosie si apare urmatorul mesaj: "No solution was found Try changing the guess value or the value of TOL or CTOL " (Solutia nu a fost gasita Schimbati valoarea initiala sau valorile TOL sau CTOL ) Fig 15 Din graficul functiei se vede ca ecuatia are doua solutii Pentru a le obtine este suficient sa schimbati valoarea initiala a variabilei Procedura given-find foloseste pentru determinarea solutiei trei metode, care sunt diferite modificari a metodei lui Newton Alegerea metodei de rezolvare se face automat Pentru verificarea solutiei obtinute, sau daca nu obtineti nici o solutie, puteti schimba metoda de rezolvare Efectuati click-dreapta pe functia find; din meniul care apare alegeti optiunea Nonlinear si din submeniu metoda dorita Cu procedura given-find poate fi rezolvat un sistem de pâna la 200 ecuatii 1 11 ECUATII ALGEBRICE LINIARE Un sistem de n ecuatii algebrice liniare cu n necunoscute are urmatoarea forma dezvoltata: (1) unde aij (i=1,n; j=1,n;) - coeficientii necunoscutelor; xj - necunoscutele sistemului; bi - termenii liberi ai sistemului Sistemul (1) poate fi prezentat si în forma matriceala: A·x=b (2) Aici A - matricea coeficientilor sistemului cu dimensiunile n×n, x - vectorul necunoscutelor cu dimensiunea n, b - vectorul termenilor liberi cu aceeasi dimensiune Matricea A se numeste singulara, daca determinantul ei este zero (|A|=0) În cursul de algebra liniara se demonstreaza, ca daca matricea A coeficientilor sistemului de ecuatii (1) nu este singulara, sistemul de ecuatii are o singura solutie În acest caz ecuatiile sistemului sunt independente, iar sistemul este compatibil Daca matricea A este singulara si vectorul termenilor liberi nu este nul (b?0), sistemul de ecuatii (1) n-are solutii si se numeste incompatibil Daca b=0, sistemul de ecuatii (1) se numeste omogen, în caz contrar neomogen Pentru rezolvarea sistemelor de ecuatii algebrice MATHCAD propune câteva instrumente: 1 procedura given-find; 2 functia lsolve 3 cu operatii matriceale directe Procedura given-find este cea mai clara, dar este comoda doar pentru sisteme cu un numar mic de ecuatii Pentru sisteme mari de ecuatii se folosesc operatii directe matriceale sau functia lsolve Pentru calculul numeric cu procedura given-find este necesar de a indica valori arbitrare initiale ale variabilelor, iar pentru calcul simbolic nu este necesar Functia lsolve are sintaxa: lsolve(A,b) Unde A - matricea coeficientilor sistemului de ecuatii, b - vectorul termenilor liberi Prezentam aplicarea ei si a operatiilor matriceale directe pentru rezolvarea sistemelor de ecuatii liniare: Daca sistemul de ecuatii algebrice este omogen si matricea coeficientilor A este singulara, sistemul are o infinitate de solutii Aceste solutii sunt infinitatea de vectori coliniari cu vectorul propriu al matricei A, care corespunde valorii proprii egala cu zero 1 12 ECUATII DIFERENTIALE Ecuatiile, ce contin functia necunoscuta sub semnul derivatei se numesc ecuatii diferentiale Deosebesc ecuatii diferentiale obisnuite si cu derivate partiale, liniare si neliniare, omogene si neomogene, cu conditii initiale si cu conditii la limita Programul MATHCAD rezolva majoritatea ecuatiilor obisnuite cu conditii initiale (asa numita problema lui Cauchy), multe ecuatii obisnuite cu conditii la limita, si chiar unele ecuatii cu derivate partiale În timpul apropiat se asteapta aparitia pe piata a unei versiuni noi a programului cu posibilitati avansate de rezolvare a ecuatiilor diferentiale Functiile, care permit rezolvarea ecuatiilor diferentiale obisnuite sunt accesibile cu comanda Insert, Function Pe ecran apare fereastra de dialog "Insert Function" Selectati din lista stânga Differential Equation Solving În lista din dreapta obtineti toate functiile disponibile cu o descriere a destinatiei, argumentilor si modului de aplicare Rezolvarea ecuatiilor se face cu aplicarea metodelor numerice Programul nu prevede (pâna ce) solutionare analitica a ecuatiilor diferentiale O singura ecuatie se rezolva cel mai simplu cu procedura Given-Odesolve, care foloseste metoda numerica Runge-Kutta de ordinul 4 În fig 16 este prezentata rezolvarea unei ecuatii liniare neomogene de ordinul doi cu aplicarea acestei proceduri Functia Odesolve are trei argumenti: primul indica variabila independenta (x); al doilea indica valoarea variabilei x pâna la care se vor efectua calculele si al treilea - numarul de pasi în care metoda Runge-Kutta va rezolva ecuatia Ultimul argument este optional (poate sa nu fie indicat) Semnul derivatelor si egalul logic se introduc cu ajutorul paletelor "Calculus" si "Boolean"(fig 3) Semnul derivatei "' " se include cu combinatia de tasre + Cu functia Odesolve se rezolva numai ecuatii cu conditii initiale! (conditii la capatul stâng al intervalului de integrare) Numarul conditiilor initiale trebuie sa fie egal cu ordinul ecuatiei diferentiale si sa contina valoarea functiei si a derivatelor ei pâna la ordinul n-1 la capatul stâng al intervalului de Fig 16 integrare Cu aceasta functie este posibil de a rezolva si ecuatii diferentiale neliniare În figura alaturata este prezentata rezolvarea ecuatiei, care descrie miscarea oscilatorie libera a pendulului matematic în cazul deplasarilor mari Aici T - este unghiul de deviere a pendulului de la pozitia verticala De obicei, rezolvarea ecuatiilor diferentiale neliniare necesita cunostinte a fenomenului, care este studiat Lipsa acestor cunostinte deseori conduce la mesajul "Could not find a solution" (Solutia nu este gasita) O metoda alternativa de rezolvare a ecuatiilor diferentiale cu conditii initiale este efectuata cu functia rkfixed, care este utila la rezolvarea sistemelor de ecuatii diferentiale obisnuite Este interesant faptul, ca ecuatia de ordinul "n" poate fi redusa la un sistem de "n" ecuatii de ordinul 1 Aplicam aceasta functie pentru rezolvarea problemei prezentata în fig 16 Aici, D(x,y) este functia, care determina sistemul de doua ecuatii y'=D(x,y) la care este redusa ecuatia din fig 16 Intervalul de integrare este între 0 si 30, y0 este vectorul conditiilor initiale M este numarul de pasi în care este divizat intervalul de integrare pentru a rezolva problema Rezultatul este o matrice cu dimensiunile (M+1)×(n+1), n este ordinul ecuatiei initiale Daca efectuati click pe matricea rezultat, aveti acces la toate elementele ei Prima coloana prezinta valorile argumentului, a doilea - valorile functiei, si a treilea - valorile primei derivate Tot aici sunt prezentate graficele functiei si a portretului fazic al ecuatiei Dupa cum am mentionat programul MATHCAD permite de a rezolva si unele ecuatii diferentiale cu conditii la limita Functiile destinate acestui calcul sunt sbval si bvalfit Ecuatii cu conditii la limita sunt acelea, pentru care o parte de conditii sunt date la capatul stâng, iar celelalte în alte puncte ale intervalului de integrare Functia sbval rezolva ecuatii, când conditiile sunt cunoscute la capetele intervalului de integrare Functia bvalfit se aplica, când conditiile sunt cunoscute în trei puncte ale intervalului de integrare (la capete si într-un punct intermediar) De fapt, aceste functii nu rezolva ecuatia cu conditii la limita, ci efectueaza conversia ei la ecuatie cu conditii initiale si apoi se rezolva ecuatia cu functia rkfixed Sintaxa functiei sbval este: sbval(v,x0,x1,D,loud,score) unde: v - vector , cu valori initiale de încercare nespecificate de la capatul stâng al intervalului de integrare; x0, x1 - punctele de început si sfârsit al intervalului de integrare; D - are aceeasi semnificatie ca si în functia rkfixed; loud - vector, cu toate conditiile initiale (de la capatul stâng ) necesare Valorile necunoscute sunt valorile de încercare ale vectorului v; score - vector, de aceeasi dimensiune cu vectorul v Fiecare componenta este diferenta între o conditie cunoscuta de la capatul din dreapta a intervalului de integrare notata prin numele ei (w0, w1, ) si valoarea ei; Functia returneaza un vector cu conditiile initiale ale problemei, care nu erau în prealabil cunoscute Exemplu: Sa se rezolve ecuatia y''''=5*x cu urmatoarele conditii la limita: y(0)=0, y''(0)=2, y(2)=3, y'(2)=4 Intervalul de integrare este de la Deci x0=0 si x1=2 La capatul stâng nu sunt cunoscute doua conditii initiale Deci, vectorul v va contine doua valori de încercare(prima aproximatie), v0 si v1 Aceste valori pot fi arbitrare sau în apropierea celor asteptate Formam functia loud(x0,v), care este un vector cu patru componente (dupa numarul de conditii initiale necesare pentru rezolvarea problemei) Fiecare componenta este conditia initiala corespunzatoare, începând cu y(x0), y'(x0), y''(x0) si y'''(x0) Deoarece y'(x0) si y'''(x0) nu sunt cunoscute, în locul lor se scrie v0 si v1 respectiv(nu valorile numerice de încercare!) Formam functia score(x1,w), în conformitate cu regulele expuse anterior, luând în considerare, ca w0, w1, w2 etc sunt respectiv valorile functiei, primei derivate, derivatei a doua, etc la capatul din dreapta al intervalului de integrare Crearea functiei D(x,y) este putin mai complicata Ecuatia diferentiala initiala de ordinul n trebuie redusa la un sistem de n ecuatii de ordinul 1 Pentru aceasta notam: y0 = y, y1 = y', y2 = y'', si y3 = y''' Alcatuim sistemul de ecuatii: y0' = y1 y1'= y2 y2'= y3 y3' = 5*x Partea dreapta a acestui sistem este vectorul, ce defineste functia D(x,y) Mai departe este simplu Definim functia sbval Afisam valoarea ei Aplicam functia rkfixed pentru obtinerea functiei y(x) si derivatelor ei Trasam graficele acestor functii pentru intervalul de integrare Cu ajutorul cunostintelor obtinute este relativ simplu de a rezolva asa numita problema valorilor proprii pentru ecuatii diferentiale obisnuite liniare omogene În Rezistenta Materialelor acestea sunt problemele calculului fortelor critice la flambaj si frecventelor proprii de vibratie ale barelor Pentru rezolvarea acestor probleme se folosesc aceleasi functii sbval si bvalfit În calitate de exemplu vom studia flambajul unei bare drepte simplu rezemate cu rigiditatea la încovoiere constanta de-a lungul barei Lungimea barei este 1m Problema se reduce la rezolvarea ecuatiei diferentiale omogene de ordinul doi: y''+?*y = 0 Aici, ? - parametrul(necunoscut), care determina forta critica; y(x) - functia deplasarilor Pentru aceasta problema conditiile la limita sunt: y(0) = 0 si y(1) = 0 Ecuatiile diferentiale liniare omogene se caracterizeaza prin faptul, ca: 1 functiile y(x) se determina numai pâna la o constanta arbitrara aditiva ne egala cu zero; 2 problema are solutii nenule numai pentru un sir infinit de valori discrete a parametrului ? Astfel, problema trebuie rezolvata cu trei conditii: doua conditii initiale si o valoare a parametrului ? Dar deoarece functia y(x) se determina doar pâna la o constanta, valoarea derivatei ei la capatul stâng al barei poate fi luata arbitrara, de exemplu egala cu unitatea Acum pentru a satisface conditia de la capatul drept vom varia parametrul ?! Notam: y0 = y, y1 = y', y2 = ? Alcatuim sistemul de trei ecuatii: y0' = y1 y1'= -y2*y0 y2' = 0 Parametrul ? este constant pentru orice x si deci derivata lui este zero Acum aplicam pentru calculul lui functia sbval, asa cum am procedat în exemplul precedent: Scriem valorile argumentului x0 si x1; Aplicam parametrului ? o valoare initiala arbitrara ?0 în apropiere de zero, pentru a obtine prima valoare proprie (Reamintim, ca problema are o infinitate de valori proprii); Definim finctiile D(x,y), load si score; Calculam cu functia sbval prima valoare proprie ? Pentru a determina functia proprie care corespunde primei valori proprii aplicam functia rkfixed Pentru a obtine alte valori proprii este suficient sa introduceti alte valori de încercare Daca luati ?0 = 50, veti obtine a doua valoare proprie 22*p2, etc Câteva exemple de rezolvare a problemelor de valori proprii le gasi-ti în compartimentul "Eigenvalues and Eigenfunction" a centrului de resurse(Resourse Center) din meniul Help Functiile aplicate anterior pentru rezolvarea ecuatiilor diferentiale folosesc metoda numerica Runge-Kutta Dar, exista ecuatii care necesita un pas extrem de mic a metodei numerice Pentru rezolvarea lor metodele standard nu pot fi aplicate Aceste ecuatii se numesc rigide (stiff) Indice a ecuatiei rigide este iacobianul functiei vectoriale D(x,y) Cu cât matricea iacobianului este mai aproape de singulara, cu atât sistemul de ecuatii este mai rigid Pentru rezolvarea acestor ecuatii se folosesc functiile Stiffb si Stiffr Descrierea argumentilor acestor functii le gasiti la chemarea lor cu comanda Function a meniului Insert 1 13 ANALIZA DATELOR EXPERIMENTALE Cel mai frecvent, un set de date experimentale sunt prezentate în forma de tablou, care consta din perechi de date (xi,yi) Problema consta în aproximarea dependentei discrete yi(xi) cu dependenta continua (functia continua) y(x) În dependenta de scopul cercetarii deosebesc trei tipuri de aproximari: 1 interpolarea si extrapolarea(predictia) datelor; 2 regresia datelor; 3 filtrarea datelor La interpolare, functia y(x) trece prin punctele (xi,yi) si aproximeaza dependenta yi(xi) doar în interiorul intervalului ce contine valorile xi La extrapolare, se aproximeaza aceasta dependenta în afara intervalului ce contine toate valorile xi La regresie, functia y(x) nu trece obligatoriu prin punctele (xi,yi) Uneori tehnica de regresie se mai numeste netezirea datelor experimentale Regresia se foloseste atunci, când legea de variatie a datelor este o functie sau o combinatie de functii calitativ cunoscute De exemplu, o functie liniara, parabolica, hiperbolica, exponentiala, trigonometrica etc La filtrarea datelor, unele date (care se considera gresite sau inutile) sau se exclud din setul initial sau se reduce influenta lor în corespundere cu un oarecare algoritm de filtrare Asta face posibil micsorarea erorilor de masurare si evidentierea datelor corecte Pentru filtrarea datelor în MATHCAD cu succes se folosesc transformarile integrale Fourier si Laplace Functiile pentru efectuarea analizei datelor si descrierea argumentilor lor le gasiti si le înserati cu comanda Function al meniului Insert Din fereastra de dialog Insert Function(vezi paj 58), în partea stânga alegeti urmatoarele compartimente: Interpolation and Prediction - pentru interpolari si extrapolari; Regression and Smoothing - pentru regresii si filtrare În partea dreapta sunt afisate functiile disponibile, iar putin mai jos descrierea argumentilor functiei selectate si forma rezultatului returnat Exemple de aplicare ale lor le gasiti în meniul Help prin specificarea comenzii Resourse Center În rezultat vor fi afisate compartimentele Centrului de resurse (Fig 16) Efectuati click pe QuickSheets and Reference Tables Din cuprins va alegeti capitolul Data Analyses si daca efectuati click pe denumirea lui obtineti acces la diferite metode de analiza a datelor Probabil ati observat, ca centrul de resurse contine o multime de exemple pentru toate problemele care pot aparea în procesul de lucru cu MATHCAD Fig 16 Aici vom prezenta doar un exemplu de interpolare liniara a datelor si unele sfaturi generale La dimensionarea barelor drepte solicitate la flambaj este necesara cunoasterea dependentei coeficientului de reducere a tensiunii admisibile de coeficientul de zveltete a barei f = f(?) Aceasta dependenta este prezentata în forma de tabel Deci dependenta este discreta fi(?i) Este mai comod de a o avea în forma continua Pentru otel 3 tabelul contine 28 perechi de date Daca definiti datele printr-o matrice, sablonul careia se va crea cu comanda Insert, Matrix, matricea va ocupa un spatiu mare pe ecran si pe hârtie( daca veti dori s-o documentati) Este posibil de a rezolva aceasta problema Procedati în felul urmator: 1 Executati comanda Insert, Component Pe ecran va fi afisata fereastra Component Wizard 2 Din lista acestei ferestre de dialog apasati pe Input Table 3 Apasati Finish Veti obtine sablonul tabloului(matricei) de introducere a datelor în forma restrânsa; cu doua rânduri si doua coloane 4 Efectuati click pe sablon Dati denumire tabloului(de exemplu date) Introduceti datele(în prima coloana valorile ?i, în a doua coloana valorile fi) Atentie! Datele argumentului (aici ?i ) trebuie introduse în ordine crescatoare Pentru ca notatiile rândurilor si coloanelor sa înceapa cu unitatea scrieti putin mai sus de tablou ORIGIN = 1(Fig 17) Trecerea de la o celula la alta a tabloului se face cu tastele (în jos) si (la dreapta) Definim prima coloana cu denumirea x (nu ?!), iar a doua coloana cu y (nu f!) x si y sunt vectori cu 28 de componente fiecare Afisam acesti vectori în forma transponata În fig 17 sunt vizibile doar ultimele componente ale lor, dar daca efectuati click pe ei aveti acces la toate elementele Fig 17 Acum efectuam cea mai simpla interpolare, interpolarea liniara cu aplicarea functiei linterp(x,y,?), asa cum este indicat în fig 17 Primul argument este coeficientul de zveltete definit cu vectorul x, al doilea este coeficientul de reducere a tensiunii admisibile definit cu vectorul y,iar al treilea defineste alta denumire a variabilei independente ? Alte denumiri a functiei si argumentului sunt necesare pentru a deosebi datele discrete de cele continue Trasam datele initiale (xi,yi) si graficul functiei obtinute f(?) Pentru a efectua alte tipuri de analiza a datelor folositi Resourse Center 1 14 PROGRAMARE MATHCAD este destinat celor, care doresc sa rezolve probleme matematice fara a cunoaste programarea Dar, pentru probleme complicate sau pentru automatizarea unor proceduri, elaboratorii de la Microsoft au creat un limbaj de programare mic, dar destul de inteligent si simplu în aplicare Operatorii acestui limbaj sunt indicati în paleta Programming (fig 3) Ei pot fi introdu-si în liniile de program numai prin apasarea butoanelor corespunzatoare a acestei palete Selectarea manuala va conduce la mesaje de eroare Liniile de cod se obtin simplu: prin apasarea butonului Add Line Stergerea liniilor se face cu comenzile Delete si Back Spase de la tastiera, în dependenta de pozitia si orientarea semnului de înserare (cornierul de culoare albastra) Una din notiunile de baza a limbajului este notiunea de variabila locala Variabila locala se defineste si i se atribuie valoare cu butonul "?" al paletei " Programming"si este "vizibila" numai în interiorul corpului programului În exemplul prezentat variabila "c" se calcula cu ajutorul unui program care consta din doua linii de cod În prima linie variabilei locale "a" i se atribuie valoarea "3" În linia a doua variabilei locale "b" i se atribuie valoarea "a+5" De obicei, rezultatul ultimei linii a programului este atribuit variabilei sau functiei care defineste programul Putin mai sus de program a fost definita variabila globala cu acelasi nume "a", dar cu valoarea "2" Putin mai jos de program am afisat valoarea variabilei "a" Deci, variabila locala "a" este "vizibila" numai în interiorul programului Acelasi lucru este valabil si pentru variabila "b" Liniile de cod sunt comode pentru definirea functiilor cu discontinuitati Aici este aplicat operatorul if de pe paleta Programming Pentru efectuarea calculelor de rezistenta a constructiilor este necesar de a determina maximul global al functiei de o singura variabila pe un interval al argumentului MATHCAD-ul are doua functii de calcul a valorilor maximale Maximize si a valorilor minimale Minimize Dar, cu regret, ele determina doar valorile extremale locale Functii predefinite MATHCAD pentru calculul valorilor globale pe un interval al functiei nu exista Vom crea un program-functie , care face acest lucru Fie denumirea acestei functii maximum Rezultatul returnat de ea va fi un numar(scalar) Functia va avea trei argumenti: functia f în forma analitica, maximul global al careia se determina, valoarea argumentului la capatul stâng a si respectiv la capatul drept b al intervalului de cercetare Exista mai multe metode de determinare a valorilor extremale a functiei Noi vom aplica metoda cea mai simpla (este posibil ca ea sa fie si cea mai sigura): metoda scanarii functiei pe tot intervalul cercetat cu un pas foarte mic Memorarea tuturor valorilor functiei într-un vector Cu functia predefinita max determinam valoarea maxima a vectorului obtinut, care si este maximul global al functiei (în sens algebric) pe intervalul (a,b) Programul este prezentat în exemplul ce urmeaza Aici se aplica operatorii for si break care au aceeasi semnificatie ca si în limbajele traditionale de programare Cu modificari mici se determina si minimul absolut al functiei Este posibil sa creati un alt program care rezolva aceasta problema cu un numar mai mic de linii de cod Alte exemple de aplicare a programarii în mediul MATHCAD sunt disponibile în Centrul de resurse Resourse Center 2 APLICATII 2 1 TRANSFORMARI LINIARE DE COORDONATE Fie un sistem de coordonate cartezian si ortogonal de dreapta cu axele x1,x2,x3 Baza ortonormata asociata lui sunt vectorii (unitari) e1,e2,e3 Matricele tensorului de ordinul unu (vectorului V) si tensorului de ordinul doi T raportate la acest sistem de coordonate sunt cunoscute Sa se determine matricele acestor tensori în raport cu sistemul de coordonate obtinut din primul prin rotirea lui cu unghiul a în jurul axei x1 Reamintim, ca pentru sistemele de coordonate de dreapta, unghiurile de rotire sunt pozitive, daca sunt orientate în sens antiorar uitându-va din directia pozitiva a axei de rotire Sa se efectueze si transformarea inversa Rezolvare: Transformarile de coordonate se fac cu ajutorul matricelor de trecere de la un sistem de coordonate la altul În cazul tridimensional matricea de trecere este o matrice patrata cu dimensiunea (3×3) Fiecare linie a ei prezinta consecutiv componentele vectorilor ortonormati e1a, e2a, e3a ai noului sistem de coordonate în raport cu cel initial Stabilim numerotarea indicilor matricelor începând cu 1 Fie Determinam vectorii de baza e1a, e2a, e3a în sistemul nou de coordonate: Cream matricea de trecere de la sistemul initial la sistemul actual de coordonate: Fie: Efectuam transformarile: în forma tensoriala în forma matriceala Aici Va si Ta sunt matricele tensorilor V si T în sistemul nou de coordonate Efectuam transformarile inverse: Mentionam ca matricele de trecere de la un sistem de coordonate la altul sunt matrice ortogonale Pentru ele este caracteristic, ca determinantul lor este egal cu 1, iar produsul scalar al oricarei linii(coloane) cu alta linie(coloana) este zero Asta este evident deoarece liniile matricei sunt întocmai vectori unitari si reciproc perpendiculari 2 2 ANALIZA STARII DE TENSIUNI ÎN JURUL UNUI PUNCT MATERIAL Starea de tensiuni în jurul unui punct al corpului este determinata de tensorul tensiunilor T, cu matricea Tij (i,j=1,3) raportata la sistemul de coordonate cartezian si ortogonal x1,x2,x3 Sa se determine: 1 invariantii tensorului tensiunilor T; 2 tensiunile principale si suprafetele principale ale tensorului tensiunilor T; 3 vectorul tensiunii t care actioneaza pe suprafata definita de vectorul normalei exterioare n =ni*ei si trece prin acest punct; 4 tensorul sferic T0 si deviatorul tensiunilor s al tensorului T 5 în ce stare se afla materialul: reversibila(elastica) sau ireversibila(plastica) în acest punct, daca materialul este otel 3 (STAS 380-71) cu limita de elasticitate sy=250 MPa Se va aplica criteriul de rezistenta von Mises (criteriul energiei potentiale specifice de variatie a formei); 6 suprafata caracteristica a tensorului T Rezolvare: Fie 2 2 4 3 74 Initiere în MATHCAD 75 Prezentarea programului 